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Abstract 

Does a given system of linear equations Ax = b have a nonnegative integer solution? This 
is a fundamental question in many areas, such as operations research, number theory, and 
statistics. In terms of optimization, this is called an integer feasibility problem. A generalized 
integer feasibility problem is to find b such that there does not exist a nonnegative integral 
solution in the system with a given A. One such problem is the well-known Frobenius 
problem. In this paper we study the generalized integer feasibility problem and also the 
multi-dimensional Frobenius problem. To study a family of systems with no nonnegative 
integer solution, wc focus on a commutative semigroup generated by a finite subset of Z'^ 
and its saturation. An clement in the difference of the semigroup and its saturation is called 
a "hole" . We show the necessary and sufficient conditions for the finitencss of the set of holes. 
Also we define fundamental holes and saturation points of a commutative semigroup. Then, 
we show the simultaneous finiteness of the set of holes, the set of non-saturation points, and 
the set of generators for saturation points. As examples we consider some three- and four- 
way contingency tables from statistics and apply our results to them. Then we will discuss 
the time complexities of our algorithms. 



Key words and phrases: contingency tables, data security, Frobenius problem, indis- 
pensable move, Markov basis, monoid, Hilbert basis, linear integer feasibility prob- 
lem, saturation, semigroup 



1 Introduction 

Consider the following system of linear equations and inequalities: 

Ax = h, > 0, (1) 

where A G Z'^^" and h £ Z'^. Suppose the solution set {x £ W : Ax = b, x > ^ %. The 
linear integer feasibility problem is to ask whether the system in ([T]) has an integral solution or 
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not. A generalized integer feasibility problem is to find all b such that there does not exist a 
nonnegative integral solution in the system with a given A. Note that there exists an integral 
solution for the system in ([1]) if and only if b is in the semigroup generated by the column vectors 
of A. Prom this, we can write this problem as follows. 

Problem 1.1. Let ai, . . . , a„ G be columns of A and 



Q = Q{A) = {aixi H h a„2;„ : xi, . . . ,Xn G Z+} 



(2) 



be the set of all nonnegative integer combinations o/ai, . . . , a„ or in other words the semigroup 
generated by ai, . . . , an- Compute a finite representation of all vectors of Q. 



Barvinok and Woods! (120031 ) introduced an algorithm to encode all vectors in the semigroup 



Q into a generating function as a short rational generating function in polynomial time when d 
and n are fixed. Therefore, using this algorithm one can compute a finite representation of all 
vectors not in Q in polynomial time if we fix d and n. However, their algorithm is yet technically 
difficult to implement so that we do not know whether it is practical or not. Modifying Problem 
II. H in this paper, we would like to solve the following problem. 

Problem 1.2. Let Q be defined as in Problem \L1[ Decide whether there is a finite number of 
integral vectors not in Q but in its saturation. 

In other words, for fixed A, decide whether there is a finite number of integral vectors b € 
such that the system in ([T]) has a nonnegative rational solution but not a nonnegative integral 
solution. 

Intensive research has been carried out on integer feasibility problems. In 1972, iKarpI (|l972l ) 
showed that solving the integer linear feasibility problem is NP hard. In the 1980's, H.W. 
Lenstra, Jr. d eveloped an algorithm, t o detect intege r solutions in the system ([T]) using the 
LLL-algorithm 



Grotschel et al. 



(|l993l ): iLenstral (|l983l )]. Lenstra also showed that integer pro- 
gramming problems with a fixed number of variables can be solved in time polynomial in the 
input size. The algorithm was actually developed in order to prove that the integer feasibility 
problem can be solved in polynomial t ime if the dimension is fixed. A late r algor ithm of similar 



Lovasz and Scarf 



structure, by 



Aardal and collaborators 



(119921') ■ was impleni e nted by 



Aardal and Lenstral (|2002l ): 



Cook et al 



Aardal et al 



(2002 



|l993|). In addition, 
2000)] have used the 



LLL-procedure to rewrite a system of linear equations into an equivalent system that was easier 
to solve with the branch-and-bound method for testing integer feasibility. In the 1990's, based 
on work by the geometers Brion, Khovanski, Lawrence, and Pukhlikov, Barvinok discovered an 
algorithm to count integer points i n rational polyto p es, and this algorithm also runs in p olyno- 
mial time if we fix the dimension [Barvinokl ()1994l ): iBarvinok and PommersheimI (jl999l )]. The 



idea of the algorithm is to encode all the integer solutions for the system in ([T]) into a rational 
generating function. 

In recent years, the generalized integer linear feasibility problem has found applications in 
many research areas, such as number theory and statistics. One such problem is the well-known 



2 



Frobenius problem, that is, for d = 1 and relatively prime positive integers oi, . . . ,a„, it is to 
ind the biggest pos i tive i nteger b such that there does not exist an integral solution in ([T|) 



Aardal and Lenstral (|2002l )]. Equivalently, it is to find the smallest positive integer b' such that 



there exists an integral solution with b = b' + b for any b G Z+ in ([T]). Since Georg Frobenius 
bcused o n this problem, it attracted substantial attention over more than a hundred years (see 



AlfonsinI (|2006l )] for a nice survey). We can generalize the Frobenius problem to the multi- 
dimensional case. Let ai,...,a„ € Z"^ such that the lattice L generated by them is H^. Let 
K = cone(ai, . . . , a„) be the cone generated by ai, . . . , a„ and l et Q in (|2l) be the semigrou p 



generated by ai, . . . , a„. Let S = {6 G Q : 6 + (If n L) C Q}. In [Miller and Sturmfeld (|2005l )]. 
a vector 5 G S" is called a saturation point in Q. We ask to find "minimal" elements of S. In the 
multi-dimensional version of the Frobenius problem, the notion of minimality can be defined in 
several ways. We present three definitions of minimality and show finiteness results of the set 
of the minimal elements of S for each definition. 

In statistics, one can find an a pplication in the data security problem of multi-way con- 

(j2003l )]. The 3- dimensional integer planar transportation problem 



Dobra et al. 



tingency tables 

(3-DIPTP) is an integer feasibility problem which asks whether there exists a three dimensional 
contingency table with the given 2-marginals or not. (In graph theory, a graph is called planar 
if it can be dra wn in a p l ane w ithout graph edges crossing.) For more details on the 3-DIPTP, 



see 



Coxl (j2002l )]. IVlachI (jl986l ) provides an excellent summary of attempts on 3-DIPTP. 



The linear integer feasibili ty problem is also closely related to the theory of Markov bases 



Diaconis and Sturmfelsl (119981 )] for sampling contingency tables with given marginals by Markov 



chain M onte Carlo methods . The notion of indispensabl e moves of Markov ba ses was de- 
fined in Takemura and Aokil (j2004l )] and further studied in Ohsugi and Hibil (j2005l )]. Recently 
Ohsugi and Hibil ( 2nnfil )" gave a simple explicit method to construct infeasible equations of ([I]) 



from non-squarefree indispensable moves of Markov bases. One finds more details in a discussion 
of three-way tables in Section [5j 

In Section [2] we define saturation points and then we will state our main theorem. Theorem 
2.51 which shows the simultaneous finiteness of the set of holes, which is the difference between 
the semigroup and its saturation, the set of non- saturation points of the semigroup, and the set 
of generators for saturation points. In Section [3l we show the necessary and sufficient condition 
for the finiteness of the set of holes. Section [J] shows a proof of Theorem 12.51 Section [5] contains 
various computational results for three- and four- way contingency tables. Section [6] will discuss 
that (1) solving Problem 11.11 (2) solving Problem 11.21 (3) computing the set of holes, and (4) 
computing the set of fundamental holes are polynomial time in fixed d and n. 



2 Notation and the main theorem 



In this section we will remind the reade r of some definitions a.nd we will se t appropriate not ation. 



We follow the notation in Chapter 7 of Miller and Sturmfelsl (|2005l )] and 



Sturmfelf 



(119961 )]. Let 
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A G Z'^^'" and let ai, . . . , a„ denote the columns of A. Let N = Z+ = {0, 1, . . .}. 

Definition 2.1. Let Q in (0j 6e the semigroup generated by ai, . . . , a„,, let K = cone(ai, . . . , a„) 
be the cone generated by ai, . . . , a„, and let L be the lattice generated by ai, . . . , a„. Then the 
semigroup Qsat = K L is called the saturation of the semigroup Q. Q C. Qsax. oind we call Q 
saturated if Q = Qsat (also this is called normal/ H = Qsat \Q is the set of holes, a ^ Q is 
called a saturation point if a + Qsat C Q. 

We assume L = 11^ without loss of generality for our theoretical developments in Sections 3 
and 4. This is for convenience in working with the Hilbert basis of K. The following is a list of 
some notations through this paper: 

K = AR" = {aixi H V anXn ■■ Xi, . . . ,Xn ^'^+} 

Qsat = K r\L = saturation oi AZ) Q 

H = Qsat \ Q = holes in Qsat 

5 = {a G Q : a + Qsa.t C Q} = saturation points of Q 

S = Q \ S = non-saturation points of Q 

We assume that there exists c € Q'^ such that c • aj > for i = 1, . . . ,n, where • is the standard 
inner p roduct. Under this a s sump tion K and Q are pointed and S is non-empty by Problem 
7.15 of [Miller and Sturmfela (|2005l )]. Qsat is partitioned as 

Qs^t = HuSuS = HuQ. 

Equivalently 

SCQC Qsat (3) 

and the differences of these two inclusions are S and H, respectively. 

If Q is saturated (equivalently H = then € S and S = Q, because Q = + QcS + Qc 
S. Therefore S = Q = Qsat in Similarly if S = Q, then G S and Qsat C Q, implying Q is 
saturated. From this consideration it follows that either S = Q = Qsax or the two inclusions in 
([3]) are simultaneously strict. 

We now consider three different notions of the minimality of saturation points, i.e., points of 
S which are minimal with respect to Q, and Qsat- We call a G S" an S'-minimal (a Q-minimal, 
a Qsat-minimal, resp.) if there exists no other b & S , b a, such that a — bGS {Q, Qsat, resp.). 
More formally a € 5 is 

a) an S'-minimal saturation point if (a -|- {—{S U {0}))) 3 = {a}, 

b) a Q-minimal saturation point if (a -|- {—Q)) Ci S = {a}, 

c) a Qsat-minimal saturation point if (a -|- (—Qsat)) 8 = {a}. 



4 



Let min(5'; S) denote the set of 5-minimal saturation points, min(S'; Q) the set of Q-minimal 
saturation points, and min(S'; Qsat) the set of Qsat-minimal saturation points. Because of the 
inclusion ([3]), it fohows that 

min(S'; Qsat) C min(5; Q) C min(S'; S). (4) 

If a € -ff , then for any b £ Q, either a — b^ Qsat or a — 6 € H. This is because if a — 5 G Qsat 
and a — b ^ H , then a — b £ Q, and hence a = b + (a — b) £ Q, which contradicts a £ H. This 
relation can be expressed as 

Qsat n{H + (-Q)) = H. 
This relation suggests the following definition. 
Definition 2.2. We call a S Qsat, a ^ 0, a fundamental hole if 

Qsat n (a + (-Q)) = {a}. 
Let Hq be the set of fundamental holes. 

Example 2.3. Consider the OTIC- difriGTisioTKil cxcLuiplG -A — (3 5 7) with L — ^sat — 
{0,1,...}, Q = {0,3,5,6,7,...}, -Q = {0,-3, -5, -6, -7,... },// = {!, 2, 4}, 5 = {5, 6, 7, . . . } 
and S = {0,3}. Among the 3 holes, 1 and 2 are fundamental. For example, 2 £ H is funda- 
mental because 

{0,1,...} n {2, -1,-3, -4, -5,...} = {2}. 

On the other hand 4 £ H is not fundamental because 

{0,l,...}n{4,l,-l,-2,-3,...} = {4,l}. 

If ^ a £ Q, then Qsat H (a + (— Q)) D {a, 0} and a is not a fundamental hole. This implies 
that a fundamental hole is a hole. For every non- fundamental hole x, there exists y £ H such 
that X — y £ Q. If y is not fundamental we can repeat this procedure. Since the procedure 
has to stop in finite number of steps, it follows that every non-fundamental hole x can be written 

as 

x = y + a, y £ Hq, a G Q, a / 0. (5) 

We also focus on a Hilbert basis of a cone K and in the next section we will show a relation 
between the set of holes H and the minimal Hilbert basis of a pointed cone K. 

Definition 2.4. We call a finite subset B C K nZ,'^ a Hilbert basis of a cone K if any integral 
point in K can be written as a nonnegative integral linear combination of elements in B. If B 
is minimal in terms of inclusion then we call it a minimal Hilbert basis of K . 

Note that there exists a Hilbert basis for any rational p olyhedral c one a nd also if a cone is 
pointed then there exists a unique minimal Hilbert basis [see lSchriiverl (jl986l ) for more details]. 
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Now we will present our main theorem of this paper and then we will present small examples 
to demonstrate the theorem. In the theorem, cone(S') denotes the set of finite nonnegative real 
combinations of elements of S and "rational polyhedral cone" is a closed cone defined by rational 
linear weak inequalities (inequalities that permit the equality case). One can find a proof of this 
theorem in Sectional 

Theorem 2.5. The following statements are equivalent. 

1. min(S'; S) is finite. 

2. cone(S') is a rational polyhedral cone. 

3. There is some s & S on every extreme ray of K. 
4- H is finite. 

5. S is finite. 




Figure 1: White circles represent non-saturation points, a triangle represents a hole, white 
squares represent ^-minimal saturation points, and black circles represent non S-minimal satu- 
ration points in the semigroup in Example 12.61 



Example 2.6. Let A be an integral matrix such that 

A = 




The set of holes H consists of only one element {{I, 2Y}. 5 = {(0,0)*}. min(5; S") = {(1, 0)*, (1,1)*, 
(1,3)*, (1,4)*}. Thus, H, S, and mm{S; S) are all finite. 
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□ □□□□□□□ 

□ □□□□□□□ 
AAAAAAAA 

ooooooooo 



Figure 2: White circles represent non-saturation points, triangles represent holes, white squares 
represent S'-minimal saturation points, and black circles represent non S-minimal saturation 
points in the semigroup in Example 12.71 



Example 2.7. Let A he an integral matrix such that 

y 2 3 4 y 

The set of holes H are the elements {(/c, 1) : A; € Z, > 1}. S = {(i,0)* : i ^ Z, i > 0}, and 
min(5;5) = {(/c,j)* : k e Z, k > 1, 2 < j < 3} U {(1,4)}. Thus, H, S, and min(S'; 5) are all 
infinite. However, min(5; Q) = {(1,2)*, (1,3)*, (1,4)*} is finite. 



3 Necessary and sufficient condition of finiteness of a set of holes 

In this section we give a necessary and sufficient condition of finiteness of the set of holes H. 
Firstly we will show the necessary and sufficient condition in terms of the set of fundamental 
holes Hq. Then we generalize the statement, such that it is stated in terms of the minimal 
Hilbert basis of K. Ezra Miller has kindly pointed out to the authors that many of our results 
can be proved more succinctly by appropriate algebraic methods. However for the sake of self- 
contained presentation we provide our own proofs and summarize his comments in Remark 13.21 
and Remark 13.61 below. 

First we show that the set of fundamental holes, Hq, is finite. 

Proposition 3.1. Hq is finite. 

Proof. Every a € Qsat can be written as 

a = aai H h Cnttn, (6) 



7 



where Cj's are nonnegative rational numbers. (If a & H, then at least one q is not integral.) If 
ci > 1, then a can be written as 

a = {(ci - [cij)ai H h c„a„} + [cijai = a + [cijai, 

and d = a — [cijtti. Therefore 

Qsat n (a + (-Q)) D {a,d} 

and a is not a fundamental hole. In this argument we can replace ci with any Cj, i > 2. This 
shows each fundamental hole has an expression ([6]), where < q < 1, i = l,...,n. However 
fundamental holes belong to a compact set. Since the lattice points in a compact set are finite, 
-f^o is finite. □ 

Remark 3.2. For any field k, consider the semigroup rings k[Q] and k[Qsa.t]- Define M = 
k[Q sat] /k[Q], which is finitely generated as a module over k[Q]. Hq is the set of degrees for the 
minimal generators of M and therefore Hq is finite. 

Let Hq = {yi, . . . , ym}- Now for each yh € Hq and each define Xhi as follows. If there 
exists some A € Z such that y^ + Xui € Q, let 

Xhi = min{A G Z | y?,, + Aa^ € Q}. (7) 

Otherwise define Xhi = oo. Note that Xhi > because yh is a hole. Then we have the following 
result: 

Theorem 3.3. H is finite if and only if Xhi < oo for all h = 1, . . . ,m and all i = 1, . . . , n. 

Proof. For one direction, assume that Xhi = oo for some h and i. Then yh + Aoj, A = 1,2,..., 
all belong to Qsat but do not belong to Q. Therefore they are holes. Hence H is infinite. 

For the other direction, assume that Xhi < oo for all /i = 1, . . . , m and all i = 1, . . . , n. By 
([5]), each hole can be written as 



n 

X 

i=l 



for some h and Aj € N, i = 1, . . . , n. Now suppose that Aj > Xhi for some i. Then 

yh + Xitti £ Q 

and 

X = yh + Xitti + ^ Xjttj G Q, 
which contradicts that a; is a hole. Therefore if a; is a hole, then Aj < Xhi for all i. Then 

n 

H C {yh + Xhitti I /i = 1, . . . , m, Q<Xhi< Xm}. 

i=l 

The right-hand side is finite. □ 
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Remark 3.4. There are several remarks to make. For each I < i < n, let 
he the semigroup spanned by aj,j 7^ i. Furthermore write 

Q[i) = ^otj + Q{i)- 

For each h and i, Xhi is finite if and only if yh € Q(i)- Since yh is a hole, actually we only need 
to check 

yh G (-Ntti) +Q(i). 

But (— Nttj) + Q(j) is another semigroup, where in A is replaced by —ai. Therefore this 
problem is a standard membership problem in a semigroup. 

Also we only need to check i such that aj is on an extreme ray. By a slight abuse of 
terminology, we simply say that aj is an extreme ray if generates an extreme ray of K. 
If there are multiple columns of A on the same extreme ray, for definiteness we choose the 
smallest one, although we can choose any one of them. Assume, without loss of generality, that 
{ai, . . . , a/;}, k < n, is the set of the extreme rays. The following corollary says that we only 
need to consider i < k. 

Corollary 3.5. H is finite if and only if Xhi < 00 for all h = 1, . . . , m and all i = 1, . . . ,k. 

Proof. The first direction is the same as above. 

For the converse direction, we show that if Xhi < 00, 1 < h < m, 1 < i < k, then Xhi < 00, 
1 < h < m, k + 1 < i < n. Now any non-extreme ray a^, i > k + I, can be written as a 
nonnegative rational combination of extreme rays: 

k 

CLi = ^ (lijCij, i> k + 1. (8) 
3=1 

Let iji > denote the l.c.m. of the denominators oi qn, . . . , qik- Then multiplying both sides by 
qi, we have 

k 

QiUi = '^{qiqij)aj, qiqij G N. 
i=i 

Also note that there is at least one qij > 0, say qij^. Consider gjttj, 2qiai, SqiUi, .... Take A G N 
such that 

XqiqijQ > A/jjp. 

Then by ^ 

k 

yh + Xqitti = yh + Xqiqij^^aj^ + ^ Xqiqijaj G Q. 

□ 
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Remark 3.6. Let k[Q], k[Qsa,t] <ind M be defined as in Remark \3.S[ The number of points in 
H is the k-vector space dimension of M . H is finite if and only if M is Artinian, which proves 
Theorem \3.3[ Let A;[Qrays] denote the monoid generated by the smallest lattice points in Q on 
the real extreme rays of Q. Then k[Q] is itself finitely generated as a module over the ^[Qrays]' 
This proves Corollary 

Another important point is that we want to state Theorem 13.31 in terms of Hilbert bases. 
Let B = {bi, . . . , b^} denote the Hilbert basis of K. As above, if bi + Aa^ € Q for some A E Z 
let 

fin = min{A € Z j 6^ + Aa^ £ Q} 
and fin = oo otherwise. Then we have the following theorem. 

Theorem 3.7. H is finite if and only if fm < oo for all I = 1, . . . ,L and all i = 1, . . . , n. 

Proof. The first direction is the same as the above proofs. 

For the converse direction, assume that fin < oo for all I = 1, . . . , L and all i = 1, . . . ,n. 
Let Uh be a fundamental hole. It can be written as a nonnegative integral combination of the 
elements of the Hilbert basis 

L 

Uh = y^^ahibi- 

1=1 

Let 

L 

A = a hi fin- 

1=1 

Then by ([5]) 

L L 

yh + Attj = ^ ahibi + ( ^ ahiflii)ai 
1=1 1=1 

L 

= y^oihi{bi + Jiuai) € Q. 
1=1 

This implies X^i < oo for all h and i. □ 

As in Corollarv 13.51 it is clear that we only need to check extreme rays among a^'s. 

Corollary 3.8. H is finite if and only if fiu < oo for all I = 1, . . . ,L and all i = 1, . . . ,k. 

Remark 3.9. In summary, determining finiteness of H is straightforward. We obtain the 
Hilbert basis B of Qsat- For each b £ B \ Q and for each extreme ai, we check 

b e (-Ntti) + Q(i). 
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Example 3.10. Let A be an integral matrix such that 



A 



( 



13 4 



1111 



) 



Then B consists of 5 elements 



B = {&i = (l,0)*,62 



(1 



1)*,63 = (1,2)*,64 



(1,3)*,&5 



(1,4)*}- 



Then we can write 63 as the following: 




-(1,0)* + 2- (1,1)* 

(1.0) * -(1,1)* + (1,3)* 

(1.1) * -(1,3)* + (1,4)* 
2 -(1,3)* -(1,4)*. 



Thus, in this case, we have fi^i = 1 for each i = 1, ... ,4 and fin = 0, where I 7^ 3 and each 



of only one element {63 = (1,2)*}. 

4 Simultaneous finiteness of holes, non-saturation points, and 
minimal saturation points 

In this section we will show the simultaneous finiteness of holes, non-saturation points, and 
(S-minimal saturation points. As in the previous section let {ai, . . . ,0^}, < n, be the set of 
the extreme rays. First, we will show the following lemmas. 

Lemma 4.1. Suppose that Q is not saturated, a ^ Q is a saturation point if and only if 
a + y ^ Q for all fundamental holes y. 

Proof. If a G Q is a saturation point, then a + y £ Q for all y S Qsat- In particular a + y £ Q 
for all fundamental holes y. 

Now suppose that a € Q is not a saturation point. Then there exists y £ Qsat such that 
a + y is a hole. This y has to be a hole, because otherwise a + y £ Q. y can be written as 
y = yh -\- b for some fundamental hole yh and b £ Q. Then a + y = a + yh + h and a + y^ has 
to be hole. Therefore we have shown that if a is not a saturation point, then a + y is a hole for 
some fundamental hole y. □ 

Lemma 4.2. Suppose that Q is not saturated. Consider any column ai of A. There exists some 
Tij € N such that niai £ S if and only if Xhi < 00 in ^ for all h = 1, . . . , m. 

Proof. This follows from Lemma 14.11 If n^ai £ S, Xhi < Ui. For the other direction take 
ni = maxhXhi. □ 



i = I, . . . , 4. Thus by Theorem \3. 7\ the number of elements in H is finite. Note that H consists 
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Now we consider the following two conditions. 
Condition 1 For each a^, there exists > such that rijaj G S. 

Condition 2 For each extreme ray ai, 1 < i < k, there exists nj > such that njaj G S. 

Proposition 4.3. Condition 1, Condition 2, and the finiteness of H are equivalent. 

Proof. Condition 1 trivially implies Condition 2. On the other hand suppose that Condition 2 
holds. Then each non-extreme ai, k < i < n, can be written as ([8]). As above let > denote 
the l.c.m. of the denominators oi qn, . . . , qik and let = x ni x • • • x nk, then n^aj G S and 
Condition 1 holds. 

Now we show the equivalence between the finiteness of H and the other two conditions. 
Using Lemma l4.2| Condition 1 is equivalent to the condition in Theorem 13.31 Also Condition 2 
is equivalent to the condition in Corollary 13. 51 □ 

Now we prove Theorem 12.51 

Proof for Theorem \2.5l [H <;=^ [2j : minf S: S) is an integral generating set of the monoid 



m 



S U jo). We then apply Theorem 1.1 (b) of [Hemmecke and Weismantell (|2006l )] or Theorem 4 



Jeroslowl (119781 1]. 

[2j <^=^ [3l : If cone(S') is not polyhedral, there must be an extreme ray e of K not in cone(5), 
since K is polyhedral. Thus, e fl 5 = 0. 

If cone(S') is polyhedral, th en it is a rational polyhedron and has a finite integral generating 
set. Thus, by Theorem 8.8 in Bertsimas and Weismantel ( 20051 )] the polyhedron cone(S') con- 



tains all lattice points from its recession cone K and {K \ cone(S')) n Z*^ is finite, which in this 
case can only happen if cone(S') = K. Thus, there is a point from S on each extreme ray of K. 

El <^=^ m : The statement El is equivalent to Condition 2. Thus, the proof follows directly by 
Proposition 14.31 

m <^=^[5l : Suppose that H is finite. Then by Condition 1, it is easy to see that S is contained 
in a compact set and hence S is finite. For the opposite implication, suppose that H is infinite. 
Since Condition 1 does not hold, there exists some i such that naj S for all n G N. Then 
{ttj, 2ai,3ai, . . .} C S and S is infinite. □ 

Now we consider the generators min(S'; Q) and we prove that min(S'; Q) is always finite. 
Then by (j4|) min(5';Qsat) is always finite as well. Note that the multi-dimensional Frobenius 
problem can be stated as computing the sets min(5; Q) and min(S'; Qsat)- 

Proposition 4.4. min(5; Q) is finite. 



Proof. Note that Q is a finitely generated monoid. Consider the algebra, k[Q] := klt"'^ , • • • j t"'" 
where k is any algebra. i c field . Then k[Q] is a finitely generated /c-algebra by Proposition 2.5 in 
Bruns and Gubeladze ( 20061 )] and therefore a Noetherian ring by a corollary of Hilbert's basis 
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theorem (Corollary 1.3 in Eisenbud (jl995l )]). Since Is ■=< : (3 £ S > is an ideal in k[Q], we 
are done. □ 



A combinatorial proof of this proposition is given in Hemmecke et al.l ()2007l )]. 



Proposition 4.5. 



min(S; Q) C min(S; Q^at) + {Hq U {0}). 



(9) 



Proof. Let a G m.m{S;Q). We want to show that a can be written as a = a + 6, where 
d G min(S'; Qsat) and b G Hq U {0}. If a itself belongs to min(S'; Qsat), then take a = a and 
fo = 0. Otherwise, if a ^ min(5; Qsat)) then by definition of Qsat-minimality there exists a' £ S 
such that ^ a — a' G Qsat- If o-' ^ min(S'; Qsat)) then we can do the same operation to a' . This 
operation has to stop in finite steps and we arrive at a G min(S'; Qsat) such that h = a — d G Qsat- 
If this fo ^ Hq, then there exists c G Q, c 7^ 0, such that b — c G Qsat- Then 

a = d + b = d+{b — c) + c, 

where a G S", fo — c G Qsat- Since S + Qsat C 5, a + (fo — c) G 5. But this contradicts 
aGmin(5;Q). □ 



5 Applications to contingency tables 

An s-way contingency table of size ni x • • • x is an array of nonnegative integers v = 
{vii,...,is)^ 1 ^ ^ nj. For < r < s, an r-marginal of v is any of the (^) possible r-way 
tables obtained by summing the entries over all but r indices. In this section we apply our 
theorem to some examples including 2x2x2x2 tables with 2-marginals and 2x2x2x2 
tables with three 2-marginals and a 3-mar g inal ( [12] [13] [14] [234]). Also we apply our theorem to 
three-way continge ncy tables from 
we used normaliz 



Bruns and Koch 



vertex representa t ion w e used CDD 



Vlach 



(|l986l )]. To compute minimal Hilbert bases of cones, 



(l200lh] and to compute e ach hyperpla ne representation and 



Fukuda 



Hemmecke et al 



(j2005l )] and Irs [Avisl ([20051)]. Also we used 4ti2 



(j2005l )] to compute matrix A for the system. 



2x2x2x2 tables 

2x2x2x2 tables with 2-marginals 

First, we would like to show some simulation results with 2x2x2x2 tables with 2-marginals, 
which can be seen as the complete graph with 4 nodes K4 and with 2 states on each node. The 



semigroup of K4 has 16 generators ai, 
that 



, ai6 in Z (without removing redundant rows) such 



1000100010001000 
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Remember that the columns of the given array are the generators of the semi g roup. All of 
these vectors are extreme rays of the cone, which we verified via cddlib [Fukudal (120051 )] . The 
Hilbert basis of the cone generated by these 16 vectors contains 17 vectors fei, . . . , bn. The first 
16 vectors are the same as aj, i.e. bi = ai, i = 1, . . . , 16. The 17-th vector bn is 

bl7 = (1 1 ... 1)* 

consisting of all I's. Thus, 617 ^ Q. Then we set the 16 systems of linear equations such that: 

Pj : bixi + b2X2 H h fciea^ie = bn 



fori = 1,2,. 
we have: 



for i / j, 



16. We solved these systems via Irs and LattE 



DeLoera et al. 



t00± ]. Then 



bn = -bi + b2 + b3 + b5 + bg + biQ, 

bn = foi - &2 + ^4 + be + &10 + bi5, 
bn = bi — bs + b^ + bj + bii + bi4, 
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= &2 


+ &3 


- &4 + &8 + ^>12 + ^>13, 


bn 


= bi 


-b5 


+ &6 + &7 + &12 + bi3, 


bn 


= b2 


+ b5 


-be + bs + bn + 614, 


bn 


= 


+ &5 


- 67 + &8 + &IO + ''IS, 


bn 


= 64 


+ be 


+ 67 - &8 + &9 + biG, 


bn 


= foi 


+ bs 


- 69 + &10 + bn + 613, 


bn 


= b2 


+ 67 


+ 69 - &10 + bi2 + 614, 


bn 


= &3 


+ be 


+ 69 - &11 + bi2 + 615, 


bn 


= 64 


+ 


+ 6in + bn — bi2 + bi6 


bn 


= &4 


+ ^5 


+ 69 - 5i3 + 614 + fol5, 


bn 


= bs 


+ ^6 


+ '>io + ''la - bu + biQ 


bn 


= b2 


+ b7 


+ bn + bi3 ~ bi5 + &16 


bn 


= bi 


+ b8 


+ bi2 + bu + 615 - &16 



Thus by Theorem 13 .7^ the number of elements in H is finite. 
2x2x2x2 tables with 2-marginals and a 3-marginal 

Now we consider 2x2x2x2 tables with three 2-marginals and a 3-marginal as the simplicial 
complex on 4 nodes [12] [13] [14] [234] and with 2 states on each node. 

After removing redundant rows (using cddlib), 2x2x2x2 tables with 2-marginals and 
a 3-marginal has the 12 x 16 matrix A. Thus the semigroup is generated by 16 vectors in 
such that: 



1 











1 











1 











1 














1 











1 











1 











1 














1 











1 











1 











1 














1 











1 











1 











1 


1 





1 

















1 





1 




















1 





1 

















1 





1 














1 





1 





1 





1 
































1 





1 





1 





1 


























1 


1 


















































1 


1 


















































1 


1 
























































1 


1 





















All of these vectors are extreme rays of the cone (verified via cddlib). The Hilbert basis of 



15 



the cone generated by these 16 vectors consists of these 16 vectors and two additional vectors 
foi7 = (1 1 1 1 1 1 1 1 1 0)*, &18 = (1 1 1 1 1 1 1 1 1 1 1)*. 



Thus, 6i7, bis Q. Then we set the system of hnear equations such that: 



bixi + b2X2 H h biQXie = hn 

xi G Z_, Xi G Z+, for i = 2, • • • , 16. 

We solved the system via Irs and CDD. We noticed that this system has no real solution (infea- 
sible). This means that 

6i7 (-Nai) + 

Thus by Theorem 13.71 the number of elements in H is infinite. 



Results on three-way tables 

Results on the saturation of 3-DIPTP are summarized in Theorem 6.4 of 



Ohsugi and Hibi 



(j200^ )]. They show that a normality (i.e., Q is saturated) or non-normality (i.e., Q is not 
saturated) of Q is not known only for the following three cases: 

5x5x3, 5x4x3, 4x4x3. 

All 2 X J X K tables are unimodular and hence saturated. This means that there is no hole in 



Q, and thus a 2 x 2 x 2 exam ple in Irving and JerrumI (|1994| )] is not a hole. All 3 x 3 x J tables 
are saturated by the result of ISullivantI (j2004l ). 
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Figure 3: An example of 3 x 4 x 6 table such that the given marginal condition is a hole of the 
semigroup. 



For 3x4x6 tables with 2-marginals, 



Vlach 



(|l986l ) showed an example which has a ta- 



ble with nonnegative real entries, but does not have a table with nonnegative integer entries. 
This example can be found in Figure [3l Actually it is a particular example of Lemma 6.1 of 
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Qhsugi and Hibil (120061 )] . lAoki and Takemural (j2003l ) presents a non-squarefree indispensable 
move z = — z~ of size 3x4x6, where 2 appears both in the positive part z"*" and the 
negative part z~ . For this z there exist two standard coordinate vectors ei,e2 such that 



u 



2ei > 0, 



262 > 0. 



In this case Lemma 6.1 of 



Qhsugi and Hibil (j2006l )] proves that h = A{u + v)/2 G Qs^t is a hole 



and this corresponds to Vlach's example. 

Using Vlach's example, one can also show that 3x4x7 tables and bigger tables have infinitely 
many holes. We take the example in Figure [3l Then we embed the table in a 3 x 4 x 7 table. 
Then we put a single arbitrary positive integer c at just one place of the seventh 3x4 slice. 
This positive integer is uniquely determined by 2-marginals of the seventh slice alone (Table 1). 
Thus for each choice of c the beginning 3x4x6 part remains to be a hole. Since c is arbitrary, 
3x4x7 table has infinite number of holes. 













sum 




c 











c 




































sum 


c 











c 



Table 1: the 7-th 3x4 slice is uniquely determined by its row and its column sums, c is an 
arbitrary positive integer. 



We can generalize this idea as follows. Let Ai denote the integer matrix corresponding to 
problem of a smaller size. Suppose that A for a larger problem can be written as a partitioned 
matrix 

^Ai 
A2 

where ^3 and A/^ are arbitrary. We consider the case that for Ai there exists a hole. Now 
consider the semigroup associated with A2. We assume that there exists infinite number of 
one-element fibers for the semigroup associated with A2. This is usually the case, because the 
fibers on the extreme ray for A2 is all one-element fibers, under the condition that A2 does not 
contain more than one extreme rays in the same direction. 
Under these assumptions consider the equation 





Ml 




tA = 





A2 




Us 


aJ 
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where ti is a hole for Ai, t2 is any of the one-element fibers for A2 and is chosen to satisfy the 
equation. Then (ti,t25^3)* is a hole for each t2- Therefore there exist infinite number of holes 
for the larger problem. 

Example 5.1. Let Ai be an integral matrix such that 

Ai = 




and let A2 = (1). From Example \2.6l H consists of only one element {ti = (1,2)*} and with A2 
we can find a family of infinite number of one- element fibers, namely ■= {c}, where c is an 
arbitrary positive integer. Let t2 = c. Then we have a matrix A such that: 



A 



Note that (ti,t2)* = (l)2,c)* is a hole for each t2 = c. Thus, since c is an arbitrary positive 
integer, there exist infinitely many holes for the semigroup generated by the columns of the matrix 
A. 
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6 Time complexity 



Barvinok and Woodsl (|2003l ) introduced an algorithm to encode all integral vectors 



In 2002, 

6 € Z'^ in Problem 11.11 as a short rational generating function in polynomial time when d and 



n are fixed (Lemma 16.31 stated below). However, in a sub-step of the algorithm they use the 
Projection Theorem which is not implementable at present. Thus, we do not know whether it 
is practical or not. From Lemma 16.31 we can show that the time complexity of computing H is 
polynomial time if we fix d and n (Corollary 16. 5p . 



One might ask the time c o mplex i ty of Problemll.2[ Usi n g the results from Barvinokl (j 19941 ) 



Barvinok and PommersheimI (jl999l ): 



Barvinok and Woods ( 



2003)], we can prove that Problem 



11.21 can be solved in polynomial time in fixed d and n (Theorem 16. ip . In order to prove the 
theorem, we will use the multivariate generating function of a set X C Z"^, f{X; x). Namely, if 
X C Ij'^, define the generating function 



fiX;x) 



sex 



X 



where denotes x^^ ■ ■ ■ x^'^ with s = (si, . . . , Sd)- If X = P fi Z'^ wit h fixed d, where P is 
a rational convex po l yhedr on, or ii X = Q with fixed d and n, then barvinok (1994) and 



Barvinok and Woodsl (j2003l ). respectively, showed that f{X;x) can be written in the form of a 



polynomial-size sum of rational function of the form: 



fiX;x) 



EX 



(10) 
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Herein, / is a finite (polynomial size) index set and all the appearing data 7j G Q and Oi, f3ij € Z'^ 
is of size polynomial. If a rational generating function f{X; x) is polynomial size in the total 
bit size of inputs, then f{X;x) is called a short rational generating function. As an example, if 
P is the one-dimensional polytope [0, N], N ^ then /(P n Z; x) = 1 + x + + • • • + x^. 



/(P n Z; x) can be represented by a short rational generating function (1 



)/(l 



Theorem 6.1. Suppose we fix d and n. There is a polynomial time algorithm in terms of the 
input size to decide whether the set of holes, H , for the semigroup, Q, generated by the columns 
of A is finite or not. 

Using the generating functions, we can show that the computation of fundamental holes for 
Q can be solved polynomial time if we fix d and n. 

Theorem 6.2. Suppose we fix d and n. Suppose Q is not saturated. The set of fundamental 
holes, Hq, can be encoded in a short rational generating function in time polynomial in terms of 
the input size. 

One notes that this algorithm outputs a generating function in the form of a short rational 
generating function. Therefore this does not return an explicit representation of Hq. However, 
if one wants to enumerate all elements in Hq, one can do the following: from the proof of 
Proposition 13. H we have Hq C P H Z,'^, where 



{x G 



: X 



^Sitti, 0<5i< 1}. 



fill 



i=l 



This shows that Hq is finite and also gives a finite procedure to enumerate i^o^ 

• Compute the Hilbert basis B of cone(ai, . . . , a„) fl L. 

• Check each z £ B whether it is a fundamental hole or not, that is, compute B Hq. 

• Generate all nonnegative integer combinations of elements in P n Hq that lie in P n Z*^ 
and check for each such z whether it is a fundamental hole or not. 



For more details, see 



Hemmecke et al 



((20071)]. 

Be fore p roofs of TheoremlG.lland Theq r eml6.2l we would like to state lemmas from 
( 2003 )] and Barvinok and Pommersheim ( IQod )]. 



Barvinok and Woods 



Lemma 6.3 ((7.3) in Barvinok and Woodsl (|2003l )]). Suppose we fix d and n. Let Q = Q{A). 
Then the generating function f{Q; x) for the semigroup Q can be computed in polynomial time 
in terms of the input size as a short rational generating function in the form of (jlOh . 



Lemma 6.4 (Theorem 4.4 in [Barvinok and PommersheimI (jl999l )]). Suppose we fix d and sup- 
pose P CW^ is a rational convex polyhedron. Then the generating function f(P fl Z'^; x) can be 
computed in polynomial time in terms of the input size as a short rational generating function 
in the form of (jlOh . 
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By Lemma 16.31 and Lemma 16.41 immediately, we have the fohowing result. 



Corollary 6.5. Suppose we fix d and n. Let Q = Q{A). Then the generating function f{H;x) 
for the set of holes, H := Qsat\Q, can be computed in polynomial time in terms of the input size 
as a short rational generating function in the form of (|10p . 



Proof. Suppose we fix d and n. By Lemma [6.31 we can compute the generating function f{Q; x) 
for the semigroup Q in polynomial time and by Lemma 16.41 we can compute the generating 
function f{Qsat',x) for the semigroup Qsat in polynomial time. The generating function f{H;x) 
fov H is f{Q sat, x) - f{Q;x). □ 



Using Corollary 16.51 we can prove Theorem 16.11 

Proof of Theorem \6.1\ Suppose we fix d and n. First, we use Corollary 16.51 to compute the 
generating function, f{H; x), for H in polynomial time in the form of (llOh . Let 



Then, we will do the following: F irst we choose / € Z'^ so tha t {l,Pij) ^ 0. We find such / in 



Barvinok and Woods 



(120031 1]. Let Z = (Ai, . . . , A^) € Z^. For 

{l,Pij) and z/j = {l,ai). Then we apply 



polynomial time by Lemma 2.5 in 
T > 0, let Xr = (exp(rAi), . . . ,exp(TArf)) and let 
the monomial substitution Xi — > expfrA,;) . We can do this mon o mial substitution in polynomial 
time by Lemma 2.5 and Theorem 2.6 in 



Barvinok and Woods 



(j2003|)]. Then 



fiH;Xr) 



^ r^exp(rz^i) 

.tr^^n-=i(i-exp(re.,)). 



Let 



h^{T) 



t'^ exp(ri^j) 



nj=i(l -expKij)) 



is a holomorphic function in a neighborhood of r = and we take the Taylor expansion around 
r = (i.e., we take the Laurent expansion around r = for /ij(T)/T'^). The coefficients of the 
kth powers, where 0</c<d— 1, of the Taylor expansion of hi are: 



1 



k I 



-^tdfc-z(^ji, • • ■,S,id) 
\i=o ' / 



where tdi{(,ii, . . . , ^j^) is a homogeneous polynomial of degree I an d which is called the Ith Todd 



Barvinok and Pommersheim 



polyno mial in di, ■ ■ ■ ,S,id (see more details in Definition 5.1 in 
tS)]). 

Now we claim that if the coefficients of negative powers of the Laurent expansion of (X^ie/ hi{T))/T'^ 
are all canceled, then H has to be finite. We prove this by contradiction. Suppose H is infinite. 
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Then, since all coefficients of negative powers in the Laurent expansion are canceled, the sum 
of the coefficients of the constant terms: 



i&i \i=o / 

must b e equal to the number of elements in H when we send r ^ ((5.2) 



fl2) 



Barvinok and Pommersheim 



(|l999l )]). Thus, the sum of the coefficients of the constant terms in ()12p must be equal to infinity. 
Since / is a finite index set, a coefficient of the constant term in the Laurent expansion of some 
rational function must be infinite. However, the Todd polynomials are polynomials in C so it is 
impossible. Thus we reach a contradiction. 

Conversely, it is obvious that if the coefficients of negative powers of the Laurent expansion 
of ^'^^ i^ot canceled, then H is infinite. 

Therefore we will have to check all coefficients of the fcth powers, where < < d — 1, of the 
Taylor expansion of each /ii(r). Since we have the polynomial size index set / and we have to 
only check d coefficients for each rational function, this computation can be done in polynomial 
time. □ 

Now we would like to discuss the intersection algorithm, which we need to encode Hq in a 
short rational generating function in polynomial time in fixed d and n. 



Lemma 6.6 (Theorem 3.6 in lBarvinok and Woodsl ( 20031 )). Let Si, S2 be finite subsets ofU^ , for 



fixed d. Let f{Si;x) and /(S'2; x) be their generating functions, given as short rational generating 
functions with at most k binomials in each denominator. Then there exist a polynomial time 
algorithm, which, given f {Si; x), computes 



f{Sir\S2;x) = 



(1 - X^'l) • • • (1 - x^'-) 

with s < 2k, where the ji are rational numbers, Ui,Vij nonzero integers, and I is a polynomial- 
size index set. 

The essential step in th e intersection algorithm is the Hadamard product [Definition 3.2 in 
Barvinok and Woodsl ( 20031 )]. Using Lemma 16.61 we can compute the union of s sets in Z"^ in 



polynomial time for fixed d and s. 
We now give a proof of Theorem 

Proof of Theorem \6.2[ Suppose Q is not saturated. Using Lemma 16.41 we compute the gen- 
erating function f{P fi Z'^;x) in polynomial time, where P is given in (jlip . Note that there 
are 2" points in P HQ, namely {x ^ Q : x = J2i=i^i^i^ ^ {^i !}}• ^o we can enumerate 
all points in P Ci Q in constant time. Let H = {P n Z,'^) \ {P n Q). Its generating function 
f{H;x) is f{P n Z'^;x) — f{P PI Q;x) and it can be computed in polynomial time. Note that 
Ho = H\{{H + {Pr\ Q)) n H) from the definition of Hq and Hq C {P n Z'^) \ (P n Q). 
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We compute the generating function for (H + (P n Q)) by the fohowing: let {PriQ) \ {0} = 
{zi, ■ ■ ■ , Z2"-i}- For each i = 1, 2, . . . , 2" — 1, let gi{x) := x^* • f{H; x) which is the generating 
function for the set Zi + H. Since 2" — 1 is a constant (we are fixing n as a constant), applying 
Lemma 16.61 we can compute the generating function for the union oi Zi + H in polynomial time. 
Since H and {P DQ) are finite we are done. □ 
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